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^_ Abstract 

Q ■ Recent transport properties on the stripe phase in Lai STsBaoi 25CUO4 by Li et al?- found 2- 

o ■ 

A . dimensional superconductivity over a wide temperature range including a Berezinski-Kosterlitz- 

^ \ Thouless transition at a temperature T=16K, with 3-dimensional superconducting (SC) ordering 



only at T=4K. These results contradict the long standing belief that the onset of superconductivity 
is suppressed by stripe ordering and suggest coexistence of stripe and SC phases. The lack of 3-D 
superconducting order above T=4K requires an antiphase ordering in the SC state to suppress 



fi I the interlayer Josephson coupling as proposed by Berg et al? . Here we use a renormalized mean 



field theory for a generalized t-J model to examine in detail the energetics of the spin and charge 
stripe ordered SC states including possible antiphase domains in the SC order. We find that the 
energies of these modulated states are very close to each other and that the anisotropy present in 
the low temperature tetragonal crystal structure favors stripe resonating valence bond states. The 
stripe antiphase SC states are found to have energies very close,but always above, the ground state 
energy which suggests additional physical effects are responsible for their stability. 

PACS numbers: 71.10.-w, 71.27.+a, 74.20.-z, 74.72.-h 



I. INTRODUCTION 

Recently Li and coworkers^ reported new results on transport properties of the stripe 
phase in Lai.875Bao.i25Cu04. They found that 2-dimensional superconducting (SC) fluctua- 
tions appear at an onset temperature T^°(=42K) which greatly exceeds the critical temper- 
ature for 3-dimensional SC order, Tc (=4K). These results contradicted the long standing 
belief that the onset of SC behavior was suppressed to very low temperatures in the presence 
of the static spin and charge density wave (SDW and CDW hereafter) or more precisely spin 
and charge stripe orderings. Li et a/.,- found strong evidence for a Berezinskii-Kosterlitz- 
Thouless transition (BKT) at Tbkt (=16K). This implies that the Josephson coupling be- 
tween the Cu02 planes strictly vanishes for T>Tc. Shortly afterwards Berg et ai— proposed 
that the strict interplanar decoupling arises because the planar superconductivity contains 
a periodic array of lines of vr-phase shift which rotate through 7r/2 up the c-axis together 
with the spin and charge stripe ordering in the low temperature tetragonal (LTT) phase. 
SDW order also appears at the same onset temperature, T^^ in zero magnetic field and 
this temperature is clearly separated from the crystallographic transition temperature T^o 
separating the low temperature orthorhombic (LTO) and LTT phases. In this material the 
LTT phase shows a superlattice ordering at all temperatures below Tco-^ Note however recent 
experiments by Fink et al. on Lai.8-xEuo.2SrxCu04"^ found different temperatures with the 
superlattice onset below the crystallographic phase transition temperature. Earlier studies 
by Lee et al. on superoxygenated La2Cu04^ found the same onset temperature for both SC 
and SDW order (T=42K). They also noted that signs of a CDW superlattice at higher tem- 
perature (T=55K) has been reported. These temperatures coincide with the values found by 
Li et al. in Lai.875Bao.i25Cu04 which suggests that Lee et al. were observing a similar stripe 
order with coexisting SDW and SC. In this case, however, the SC order is 3-dimensional, 
consistent with the absence of 7r/2-rotations in the crystal structure. These experiments lead 
us to conclude that in the presence of a CDW superlattice, coexisting SDW and antiphase 
d-wave SC can be favored. 

Actually a similar ordering was suggested on general grounds earlier by Zhang^ and 
also by Himeda , Kato and Ogata^ on the basis of variational Monte Carlo calculations 
(VMC) for the strongly correlated one band t — t' — J model. Himeda et af found that 
a modulated state with combined SDW, and CDW and d-wave superconductivity (dSC) 



containing site- or bond- centered anti-phase domain walls (vrDW) ( a state we denote as 
SDW+CDW+APdSC^/'') had a lower energy than a uniform d-wave SC state over a wide 
range of parameters and was even lower than a modulated state without anti-phase (denoted 
as SDW+CDW+dSC*/^) in a narrower parameter range. Recent VMC and renormalized 
mean field theory (RMFT) calculations^ have found that CDW+APdSC^/^ state (tt-DRVB 
state in ref.-) cost surprisingly little energy even in the absence of SDW modulations. 

In this paper we report on calculations using the RMFT method to examine in greater 
detail the energetics of these novel modulated states within the generalized t — t' — t" — J 
model. This method approximates the strong correlation condition of no double occupancy 
by Gutzwiller renormalization factors and generally agrees well with full VMC calculations 
which treat the strong correlation condition exactly. The static stripe phase appears in the 
LTT phase of Lai.875Bao.i25Cu04. This crystallographic phase is entered at a temperature 
Tco (=52K >T l^) and displays a complex crystal structure which has not been fully de- 
termined to the best of our knowledge. Note that although the overall crystal structure is 
tetragonal the individual Cu02 planes do not have square symmetry. Along one (x-) axis 
the Cu-O-Cu bonds are straight but in the perpendicular direction they are buckled^. Since 
the Cu-Cu distance is required to be the same in both directions there is a compressive stress 
along the x-axis which may well be the origin of the CDW superlattice that appears at the 
crystallographic phase transition into the LTT phase. At present the detailed displacements 
inside the supercell have not been refined. In our calculations we introduce a site depen- 
dent potential shift to mimic this effect. In addition we examine the effect of the hopping 
anisotropy between x- and y-axes which results from the different Cu-O-Cu bonding in the 
X and y directions. Such anisotropy was also considered by Capello et al}^ in their work on 
stripes made from anti-phase shifts in the superconductivity. 

II. RENORMALIZED MEAN FIELD THEORY FOR THE EXTENDED t - J 
MODEL 

The t — J model was introduced in the early days of cuprate research by Anderson and 
by Zhang and Rice to describe lightly hole doped Cu02 planes^-. In this single band model 
configurations with doubly occupied sites are strictly forbidden due to the strong onsite 
Coulomb repulsion. The Hamiltonian takes the form, suppressing the constraint 



Htj = - ^ t(ij) \c[„Cj^„ + /i.c.j + ^ J(,j)S, • Sj 

+ 5^^^ni. (1) 

i 

In the first term we include hopping processes between nearest neighboring (nn) sites (de- 
noted by {i,j)), next neighboring sites (nnn) and 3rd neighboring sites (nnnn) on a square 
lattice with matrix elements t , t' , t" respectively. We will measure all energies in unit of 
to (300 meV) — a standard value for the nn hopping matrix element t. The superexchange 
spin-spin interaction between nn sites J = 0.3, and a the spin index takes the value ±. In 
addition we introduce a potential shift Vi which varies from site to site within the super- 
cell to mimic the effect of the crystallographic superlattice in the LTT crystal structure. 
The strong coupling constraint of no double occupancy is very difficult to treat analyti- 
cally. Zhang and coworkers introduced Gutzwiller renormalization factors to approximate 
the constraints^. This approximation has been shown to be quite accurate for mean field 
theories when compared to numerical evaluations by VMC of expectation values of the corre- 
sponding mean field wavef unctions, |\E'), which are exactly projected down to the constrained 
Hilbert spaced. Later the case of AF ordering was considered by Himeda and Ogata, who 
showed that an anisotropic spin renormalization term is required to reproduce the VMC 
results^. The resulting renormalized Hamiltonian is 

J2 9li,j)Aid) (4,aCj,a + h.C.^ 



H = - 






+ Y,Vini- (2) 

i 

The renormalization factors g^, g^'^^ and g'^''"^ used to evaluate a projected mean field wave- 
function depend on the local values of the magnetic and pairing order parameters and the 
local kinetic energy and hole density which are defined as follows 

rrii = (^ol^fl^o); 

X{mv = (^o|4A^I^o); 

6i = l-(^olri.l^o), (3) 



where |\l/o) is the unprojected wavefunction. The two pairing amphtudes A^j ,,)a.=± are 
treated independently to incorporate a possible triplet component. The explicit renormal- 
ization factors introduced first by Himeda and Ogata are quite complex,— and we use here 
a simpler form as follows,— 



9{i,j),(7 9i,a9j^ai 



J , 26i (1 -6i) l + 6i + alnii 



9i,a 
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(4) 

where A{i,j) = Xlo- A(«,i>,o-/2, X{i,j) = So-X(i,i>,cr/2. Since the g-factors depends on the order 
parameters, the approach by direct diagonalization of the mean field Hartree-Fock Hamilto- 
nian obtained from the Hamiltonian Eq[2] will not give the best energy of the Hamiltonian 

Et =(^o|^|^o) 

= - Yl 9li,j),J(t,j) [X{i,j),a + h.c] 



n ' n ' " A* 






(s,z s,xy ^ \ 

+ ^ l^iTii + ^ ^ J"^.^ J(i j)m,mj- (5) 

Instead, we minimize the energy with respect to the unprojected wave function |\E^o) under 
the constraints ^^ rij = A'"e, (\E'o|^o) = 1? ^e is the total electron density. That is equivalent 



to miniinizing the function 



W = (^ol H l^o) - A ((^ol^o) - 1) - /i f 5^ 



hi - N^ 



which results in the following variational relation 

5W 
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(6) 



(7) 



For an operator O with the expectation value O = (\E'o|0|^o), (^ (^o| O |\E'o) /^ (^o| 
O l^o)- Thus one obtains the following mean field Hamiltonian, 



H 
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(8) 



which satisfies the Schrodinger equation Hmf |^o) = A l^&o)- The coefficients of Hmf are 
given as 

dW 



dX{z,j),. 



dW 



dW 
driia 



-5, 



s,z s,xy ^i 



~9{i,j),J{i,j) + 
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(^ - Vi) + 7;(^Y1 3'(',j)J{i,j)m] + 



dW 

dni„ 



(9) 



with dW/dxl,^^-,^, = [^W/^Xi^,,),.]\ dW/d^l^^^^^^ = [dW/dA(i,),^]\ 5(,,,),(,,,) = 1 only when 
i and j are nn, otherwise it equals 0, the partial derivative terms [f^] in the above equations 
refer to the derivative of W with respect to the mean field O via the Gutzwiller g-factors (see 
Eq[l]). This mean field Hamiltonian Hmf in Eq[S] is then solved self-consistently. In the 
numerical calculations, we always diagonalize Hmf for a sample consisting of 257 supercells 
along the direction with periodic boundary condition unless stated explicitly otherwise. 

III. SIMPLIFIED MODEL: SITE-CENTERED ANTI-PHASE DOMAIN WALLS 
WITH D-WAVE SUPERCONDUCTIVITY 

We begin the discussion of the results with the simplest case namely site-centered anti- 
phase domain walls in a d-wave superconductor (APdSC). To this end we restrict the 
Hamiltonian to the two terms without SDW order, and solve it self-consistently without 
considering explicitly the doping dependence of g-factors, 

-Y.9'J^KAAA + h-^- (10) 

(id) 

Note that Xi,j 7^ Xi,ji ^ind Aj .,■ ^ Aij but has the same symmetry as Aj j. To keep the 
model simple, we set Xi,j = Xp independent of {i,j), and (?* = 26/ (1 + 6), g^ = 4/ (1 + 5) 
where 5 is the average doping away from half-filling. We consider first an isolated ttDW 
which lies in the center (i^ = 28) of a finite sample with open boundary condition along x 



direction and width L^ = 55. To this end we set Ajj = Ap except for the bonds along 
the domain wall which are set to zero, i.e Aj ,,|j^=j^=28 = 0. The vr-phase shift requires that 
for the two bonds {i,j) and {i',j') which are located symmetrically on the two sides of the 
domain wall, Aij\i^j^<28 = ~^i',i'liiJi>28- The change of sign at the domain wall causes 
an Andreev bound state (ABS) to appear at the chemical potential which we take as the 
energy zero. This shows up clearly when we calculate the local density of states (LDOS) 
as illustrated in Fig[T^,b]. For the case of weak coupling in Fig[T^] a clear peak appears 
in the LDOS at zero energy for sites at or near the domain wall [i^. = 27, 28, 29) while 
far away sites show peaks at the bulk gap edges and reduced values at zero energy. This 
behavior shows up also very clearly in the spatial dependence of the quasiparticle spectral 
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FIG. 1: (Color online) Local density of states (LDOS) for a simplified model (Hamiltonian Hg 
defined in Eq[10]) for an isolated site-centered anti-phase domain wall in a d-wave SC. Periodic 
boundary condition along y direction and open boundary condition along x direction are imposed. 
The width of the system along x direction is L^ = 55 with the domain wall located at site 28. 
The average doping concentration is fixed at 5 = 0.25, and Xp = 0.20. Panels (a) and (b) are for 
Ap = 0.02, and 0.08 respectively. The 14th site is halfway between the domain wall (site 28) and 
the edge, and the two sites (27, 29) are neighbors of the domain wall. A broadening factor 0.004io 
is used. 

weight. This is illustrated in Fig[2^,b] for the case of a weak and a moderate gap value of 
the pairing amplitude Ap = 0.02(0.08). The spectral weight is concentrated close to the 
vrDW at quasiparticle energies Ek — 0, but away from the vrDW for values of Ek near the 
bulk gap energy E^ = 2g^JoAp. The total energy differences between the states with and 
without vrDW for the two Ap are 0.0066to and 0.0365to, respectively. The energy cost of 
the domain wall is substantial, consistent with the creation of a LDOS peak in the center 
of energy gap. Note that for the case of a moderate gap value of Ap, the peak of LDOS 
near Ek — shows structures consistent with the development of a one- dimensional band 
of Andreev bound states which propagate along the domain wall. This can be also seen in 
the quasiparticle dispersion which is a function only of ky. 

Turning our attention to a periodic array of parallel ttDW, we focus on the case of 
period Lx = 4, relevant to the cuprates, illustrated in Fig[3]. In this case the Andreev 
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bound states on neighboring domain walls will overlap strongly leading to a more complex 
dispersion relation for the associated quasiparticle states. Note the d-wave form of the bulk 
superconductivity leads to gapless excitations in the nodal directions which in turn leads to 
stronger overlap for near nodal quasiparticles. To illustrate this more complex behavior we 
focus on a particular model which can be solved analytically. To this end we set 5 = (i.e. 
half-filling), (7* = and set Xp = ^p and g'^JXp = 1- In this case the quasiparticle dispersion 
is obtained by diagonalizing the Hamiltonian 



Hk 



A, 



Bk 



^^ ^k ^k 


/ 




^ 2 cos ky e*''^ 





^-ikx \ 


e"*''^ 2 cos ky 


^ikx 





e-'''- 


2 cos ky 


p'^'^x 


\ e'^- 


Q-ikx 


2 cos ky J 


/ e'^- 





_^-ik, \ 


f,-ik^ -2 COS ky 


p^fcx 





6-^^=- 





_eifcx 


\ -6'^=^ 


p '^f^x 


2 cos ky J 



(11) 



where X|, 



^t 



Cjf^pCi-k,ij with / = 1,2,3,4 denoting the sites inside a supercell. The 
quasiparticle dispersion takes a simple form. 



Ek = ±\ 6 cos2 ky + A± 2\ (2 + cos^ ky) - 4 sin^ 2k^. 



(12) 



For a wavevector {k^, ky) close to (vr/2, tt/2), the two quasiparticle bands close to the Fermi 
level have an anisotropic nodal structure with 



Eu 



±2JKl + 2Kl 



(13) 



where {K^, Ky) = {n/2 — kx, 7r/2 — ky). This nodal structure completely suppresses the the 
density of states (DOS) at zero energy as shown in Fig[3], and pushes the peaks in the DOS 
of the Andreev bound states away from the chemical potential. 



IV. COEXISTING ANTI-PHASE SUPERCONDUCTIVITY AND SPIN AND 
CHARGE DENSITY WAVES 

Anti-phase domain walls in a superconductor usually cost a substantial energy. The 
key question raised by the recent experimental results of Li et al^. on the static stripe 
phase is whether SDW and CDW coexisting with vrDW lead to a state with a net energy 
gain. The VMC calculations of Himeda et air found a small energy gain for a longer 
superlattice with a larger separation between vrDW within a restricted parameter range. 
Recent calculations for a 8-superlattice without SDW order by Raczkowski et ai— did not 
yield an energy gain but the energy cost to introduce ttDW was quite small. These results 
motivated us to examine a wider parameter range within a RMFT approach and look for a 
possible net energy gain in an 8-superlattice (with site-centered anti-phase domain walls) at 
a hole concentration 6 = 1/8 when coexisting SDW order and vrDW are included. A longer 
10-superlattice (with bond-centered anti-phase domain walls) state gives similar results. 
In view of the orthorhombic nature of the individual Cu02-planes in the LTT phase, we 
allowed for anisotropy in the hopping tx(y) and exchange coupling Jx(y)- Below we keep 
the nn hopping in the y-direction fixed, ty = to, and scale Jx/Jy = tl/ty. In addition the 
presence of a crystallographic superlattice in the LTT phase motivated us to examine also 
the effect of the lattice inhomogeneity by including a site dependent potential modulation, 

A. Site-centered anti-phase dSC 

The RMFT approximation yields a series of coupled nonlinear equations. An iteration 
method is used to obtain optimal values of the four order parameters: the pairing and 
hopping amplitudes, sublattice magnetization and hole density. When the solution iterates 
to a stable set of values we can conclude that a local energy minimum exists, but on occasion 
no stable solution can be found, which indicates that no local minimum exists with this 
symmetry. In general we find stable solutions for the case of coexisting CDW and SDW 
with or without vrDW. Typical patterns for an 8-superlattice are illustrated in Fig[4] with or 
without site-centered vrDW in systems where the modulation of the pairing amplitude is site 
centered. The antiferromagnetic domain wall (AFDW) coincides with the maximum hole 
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density while the vrDW appears at the minimum hole density. (In the case without SDW, 
the ttDW appears at the maximum hole density^.) 

In table [H the results for the ground state energy and local values of the order parameters 
are presented. The upper lines are for the case of nn hopping only {f = 0), with and without, 
anisotropic component in the nn hopping tj.(y). In this case t^j. = ty the results show that 
the uniform AFM+dSC state is lowest. When AFDW and the associated modulation of 
the hole density are included the resulting state (denoted by SDW+CDW+dSC"*) has an 
energy that is slightly higher. Introducing vrDW to create antiphase superconducting order 
(SDW+CDW+APdSC) raises the energy further. Anisotropy in the nn hopping narrows the 
energy differences but does not change the relative ordering of the states with and without 
ttDW. When a weak nnn hopping is added, the SDW+CDW+dSC state gains in energy 
and when anisotropy is also added this state has the lowest energy. In this case when we 
consider anisotropic nn hopping, the energy cost of introducing vrDW in the superconducting 
is further reduced to small but still positive values. A further increase in the nnn hopping 
term (shown in Fig[5^]) however does not lead to an energy gain for ttDW. The energy cost 
of vrDW remains very small but positive. 

The presence of substantial local modulations in the hole density in these states led 
us to investigate the effect of introducing a site dependent potential shift. Such a shift can 
result from the crystallographic superlattice modulation that appears at the crystallographic 
transition into the LTT state. The results in Fig[5b] show that this potential shift reduces 
the energy cost of the site-centered anti-phase domain wall and enhances the charge and 
spin modulation but still does not lead to a net energy gain for the SDW+CDW-|-APdSC^ 
state even in the most favorable case of anisotropic nn hopping and substantial nnn hopping. 
Within the RMFT the vrDW always demands an energy cost even though it may be only a 
very small amount. Bond-centered vrDW with anisotropic nn hopping and longer periodicity 
La; = 10 shows that the energy difference between these two states, with and without vrDW, 
can be also very close. 
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Et 


Bkin 
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imax 
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^n-max 


Am„x 


Ami„ 


Xmax 
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AFM+dSC 


-0.4878 


-0.3287 


-0.1593 


0.12500 


0.12500 


0.08524 


0.1142 


0.1142 


0.1903 


0.1903 


tx 


= 1.00 


dSC 


-0.4863 


-0.3428 


-0.1435 


0.12500 


0.12500 





0.1152 


0.1152 


0.1928 


0.1928 


t' 


= 0.0 


SDW+CDW+dSC= 


-0.4865 


-0.3373 


-0.1492 


0.1372 


0.1134 


0.08412 


0.1214 


0.09917 


0.2215 


0.1821 






SDW+CDW+APdSC= 


-0.4782 


-0.3292 


-0.1490 


0.1498 


0.09604 


0.1418 


0.1114 





0.2639 


0.1111 






AFM+dSC 


-0.4536 


-0.3117 


-0.1419 


0.12500 


0.12500 


0.07432 


0.08399 


0.07724 


0.2652 


0.1098 


tx 


= 0.85 


dSC 


-0.4526 


-0.3225 


-0.1301 


0.12500 


0.12500 





0.08409 


0.07593 


0.2675 


0.1122 


t' 


= 0.0 


SDW+CDW+dSC= 


-0.4560 


-0.3050 


-0.1510 


0.1737 


0.07130 


0.1720 


0.06538 


0.05154 


0.2756 


0.07822 






SDW+CDW+APdSC= 


-0.4554 


-0.3036 


-0.1518 


0.1815 


0.05752 


0.1871 


0.04479 





0.2866 


0.06392 






SDW+CDW+dSC' 


-0.4567 


-0.3002 


-0.1564 


0.1911 


0.06029 


0.1831 


0.06692 


0.05151 


0.2769 


0.05869 






SDW+CDW+APdSC'' 


-0.4563 


-0.2978 


-0.1586 


0.1988 


0.04986 


0.1941 


0.06055 





0.2841 


0.03663 






AFM+dSC 


-0.4841 


-0.3219 


-0.1622 


0.12500 


0.12500 


0.09356 


0.1149 


0.1149 


0.1893 


0.1893 


tx 


= 1.00 


dSC 


-0.4817 


-0.3372 


-0.1445 


0.12500 


0.12500 





0.1179 


0.1179 


0.1920 


0.1920 


tx 


= -0.1 


SDW+CDW+dSC= 


-0.4829 


-0.3268 


-0.1561 


0.1525 


0.1008 


0.1268 


0.1304 


0.09917 


0.2215 


0.1654 






SDW+CDW+APdSC° 


-0.4759 


-0.3202 


-0.1557 


0.1650 


0.08549 


0.1731 


0.09773 





0.2642 


0.1034 






AFM+dSC 


-0.4507 


-0.3056 


-0.1451 


0.12500 


0.12500 


0.08474 


0.07982 


0.07170 


0.2688 


0.1031 


tx 


= 0.85 


dSC 


-0.4490 


-0.3182 


-0.1308 


0.12500 


0.12500 





0.08140 


0.07155 


0.2721 


0.1054 


t' 


= -0.1 


SDW+CDW+dSC= 


-0.4539 


-0.3024 


-0.1515 


0.1750 


0.07336 


0.1751 


0.06401 


0.04880 


0.2742 


0.08047 






SDW+CDW+APdSC= 


-0.4533 


-0.3021 


-0.1512 


0.1801 


0.06318 


0.1867 


0.04286 





0.1840 


0.06953 






SDW+CDW+dSC' 


-0.4538 


-0.2991 


-0.1547 


0.1833 


0.07206 


0.1740 


0.07046 


0.05567 


0.2728 


0.07248 






SDW+CDW+APdSC'' 


-0.4533 


-0.2972 


-0.1560 


0.1890 


0.06202 


0.1860 


0.05946 





0.2810 


0.04685 






AFM+dSC 


-0.4813 


-0.3151 


-0.1662 


0.12500 


0.12500 


0.1188 


0.1086 


0.1086 


0.1866 


0.1866 


tx 


= 1.00 


dSC 


-0.4750 


-0.3303 


-0.1446 


0.12500 


0.12500 





0.1216 


0.1216 


0.1899 


0.1899 


t' 


= -0.3 


SDW+CDW+dSC= 


-0.4814 


-0.3213 


-0.1602 


0.1709 


0.09043 


0.1746 


0.1263 


0.07922 


0.2351 


0.1280 






SDW+CDW+APdSC° 


-0.4760 


-0.3236 


-0.1523 


0.1700 


0.09028 


0.2002 


0.07064 





0.2431 


0.1266 






AFM+dSC 


-0.4491 


-0.2986 


-0.1506 


0.12500 


0.12500 


0.1127 


0.06892 


0.05955 


0.2683 


0.09673 


tx 


= 0.85 


dSC 


-0.4436 


-0.3122 


-0.1314 


0.12500 


0.12500 





0.08064 


0.06819 


0.2762 


0.09695 


t' 


= -0.3 


SDW+CDW+dSC= 


-0.4523 


-0.3008 


-0.1515 


0.1774 


0.08221 


0.1883 


0.06455 


0.04278 


0.2681 


0.08799 






SDW+CDW+APdSC 


-0.4518 


-0.3017 


-0.1501 


0.1787 


0.08034 


0.1822 


0.03503 





0.2768 


0.07811 






SDW+CDW+dSC' 


-0.4513 


-0.2985 


-0.1528 


0.1789 


0.09117 


0.1638 


0.06976 


0.04589 


0.2680 


0.08756 






SDW+CDW+APdSC'' 


-0.4506 


-0.2984 


-0.1523 


0.1806 


0.08737 


0.1700 


0.04898 





0.2762 


0.06648 



TABLE I: Key results for various states obtained by selfconsistently solving the Hamiltonian in 
Eq[8] with an average hole density of 1/8. Listed are the mean field energy Et, kinetic energy E^in, 
spin-spin superexchange energy Ej and the modulation of hole concentration {6max and 5min), 
the largest antiferromagnetic moment {\m\max), pairing amplitude Amax and Amm, and Xmax and 
Xmin (^ = J2a^^^ X = Server) for various statcs including homogeneous AFM-|-dSC and dSC 
states, SDW+CDW+dSC" (L^. = 8), SDW+CDW+APdSC" (L^. = 8) and SDW+CDW+dSC* 
(L^ = 5), SDW+CDW+APdSC^ (L^ = 10) states [APdSC (dSC) stands for d-wave SC state with 
(without) anti-phase domain walls site-centered (") or bond-centered ( )]. The energy unit is Iq = 
SOOmeV, V = 0, Jx/Jy = t'^/ty, ty = ^: Jy = 0-3- Anisotropic nn hopping tends to energetically 
favor the homogeneous AFM+dSC state compared to the homogeneous dSC state. However, it 
also causes some inhomogeneous state to be energetically more favored, here for instance, the 
SDW+CDW+dSC" state. Note that to introduce anti-phase domain wall in the pairing order 
parameter in the renormalized mean field theory for the t — J model always cost energy, although 
it can be very small. 
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B. Bond-centered anti-phase dSC 

Alternative bond-centered anti-phase modulations of the pairing amplitude were consid- 
ered by several groups.-i^»i^ In the case of the 8-superlattice we did not find any stable bond 
centered solution with nonzero SDW in the doping regime around 1/8 when requiring there 
is antiferromagnetic domain wall {mi = 0). But for longer periodicity Lx=10 we found a 
stable solution. In Fig[6] a typical pattern for this long 10-superlattice with and without 
the bond-centered vrDW is illustrated. The energy cost of the APdSC'' is also positive for 
the bond-centered case but is even smaller compared with the site-centered case (see table 
H]) at some cases. 

V. SPECTRAL PROPERTIES OF THE MODULATED PHASES 

Next we examine the density of states in the modulated phases which gives us insight into 
the interplay between the SDW and SC with either dSC or APdSC order in the stripe phases. 
We restrict our considerations to the case of site-centered pairing modulation relevant for 
8-superlattice. It is instructive to calculate several density of states, starting with the local 
density of states (LDOS) 

Aj{u) = ~-y2lmGj,,{u), (14) 

a 

where Gj^a (uj) is the Fourier transform of the time dependent onsite Green's function 
Gi^a (t) = —i (TtCi^fj (t) Cj^^ (0) V The averaging of the LDOS over all sites gives 

A(a;) = l/iV,^A,(a;), (15) 

where Nc is the size of a supercell. Also of interest is the quasiparticle (QP) density of states 

N{u) = -^ [dkJ2H^-Ei), (16) 

where / denotes all the quasiparticle bands in the reduced Brillouin zone (RBZ), Vrbz is 
the volume of RBZ, k cRBZ. This latter is the density of states which determines the sum 
of the quasiparticle eigenvalues which enters the ground state energy in mean field theory. 
The results for these DOS in the various modulated phases are presented below. First we 
consider the cases of a dSC with array of vrDW and of a SDW separately and then the results 
when both orders coexist. 
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(a) anti-phase dSC 

We start with the DOS for an array of ttDW with a superlattice periodicity of 8 and an 
average hole density of 1/8. The LDOS is shown in Fig[7], for the 3 independent lattice 
sites, site 1 at the vrDW, site 3 halfway between the ttDW and the remaining equivalent sites 
2, 4. In the energy region near zero, the prominent features are a finite LDOS at all sites, 
which is largest at the center of a vrDW (site 1) and two sharp peaks (labeled as A and B) 
symmetrically placed at positive and negative energies. The finite LDOS at ii^ = implies a 
finite quasiparticle Fermi surface in this APdSC'^ state. The quasiparticle energy dispersion 
is quite complex and is illustrated in Fig [7b]. Along the high symmetry line, k^ = 0, in RBZ 
there are 3 nodal points. These expand into nodal lines for finite k^. to create two closed 
Fermi loops shown in Fig[7b.]. The two sharp peaks labeled A and B in the DOS, A {u), can 
be shown to originate from the almost fiat bands displaced away from zero energy in Fig[7l3]. 
The LDOS that appears in Fig[71i] shows clearly an enhanced DOS near zero energy which 
implies a substantial energy cost to introduce the vrDW into a uniform dSC state. 

(b) SDW 

The second case we considered is a simple SDW state in which an array of AFDW is 
introduced to create a 8-superlattice. Again the LDOS (see Fig[8]) shows finite values at 
zero energy with the largest value at the center of the AFDW (rrij = 0). As a consequence this 
SDW state is metallic. Note a uniform state would also be metallic at this hole concentration 
of 5 = 1/8. It is however very relevant that the SDW superlattice does not truncate the 
Fermi surface completely to give an insulating state, since then coexistence with d-wave 
pairing would be disfavored. Further any coexisting state would not be superconducting. 
The Fermi surface shown in Fig[8bi] consists of standing waves along ky i.e. perpendicular 
to the AFDW and two one-dimensional bands propagating along AFDW. 

(c) Coexisting SDW, CDW and dSC or anti-phase dSC 

We examine the coexisting state to look for possible synergy between the SDW and dSC 
and also to compare the two possibilities for the superconducting uniform dSC and the 
APdSC*, i.e. superconductivity without and with an array of vrDW. The favorable choice 
of the relative position of the two domain walls is to stagger the vrDW and AFDW as shown 
in Fig[l]. From Fig[9bi,b] one sees that in both cases the LDOS develops a strong minimum 
around zero energy with even a drop to zero in a narrow range around zero energy. The site 
dependence of the LDOS is weaker than in the previous cases. This strong energy minimum 
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indicates a certain synergy between the SDW and dSC which can lower the energy through 
a truncation of the finite Fermi surface that exists in the both cases separately, SDW and 
APdSC*. The energy difference in the LDOS between the two cases with and without 
ttDW, is small. But when one compares the total ground state energy, a finite energy cost 
to introduce vrDW into the superconductivity always appears. 

The strong minimum in the DOS at the Fermi level in the SDW+CDW+APdSC'^ state 
is consistent with the spectra obtained in angle resolved photoemission (ARPES) and scan- 
ning tunnelling (STM) experiments on Lai §75 Bao.i25Cu04 reported by Valla et alr^. Our 
calculations give a complex quasiparticle dispersion associated with the 8-fold superlattice 
which does not seem to be resolved in the ARPES spectra. So a more detailed comparison 
with experiment is not possible at this time but the main feature of the experimental DOS 
is reproduced in our calculations. 

VI. DISCUSSION 

Anti-phase domain wall or ttDW generally cost considerable energy in a superconductor 
because they generate an Andreev bound state at the Fermi energy due to the interference 
between refiected electrons and holes. This effect is illustrated in Fig[T^] which shows a 
peak in the LDOS centered on an isolated vrDW. In an array of parallel vrDW this DOS 
peak broadens into a 2-dimensional band due to both the propagation of the ABS along the 
vrDW, as illustrated in Fig[T]D], and the overlap of the ABS on neighboring vrDW. This leads 
to structure which can lead to a pronounced minimum in the LDOS in certain cases such 
as the case of a closely spaced array of vrDW shown in Fig[3)D]. This structure in the LDOS 
lowers the energy cost to introduce vrDW in the dSC, but leaves it still positive. For the 
period 8 supercell the modification of the DOS is less important. As illustrated in Fig[7b] 
the APdSC* bandstructure is quite complex and displays a finite Fermi surface (see Fig[7ti]). 
The resulting LDOS has a finite value at the Fermi energy which is largest at the center of 
the ttDW. 

In the case of coexisting SDW and CDW one must first consider how the effect of these 
superlattices alone. The results are presented in Fig[8] which shows a metallic state with a 
finite DOS and Fermi surface. This is important since if the SDW resulted in an insulating 
groundstate, the addition of Cooper pairing would be less energetically favorable and would 
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not change the state from insulating to superconducting. The bandstructure consists of 
standing waves in the direction perpendicular to the AFDW which are propagating in the 
direction parallel to the AFDW. 

Coexisting SDW and dSC leads to a substantial interplay between the two broken symme- 
tries. Recently Agterberg and Tsunetsugu showed that there can be a synergy between the 
two orders due to the presence of cross terms involving both order parameters in a Landau 
expansioni^. The cross term depends crucially on the relative orientation of the wavevector 
of the SDW and APdSC For the case of parallel q- vectors under consideration here (eg. as 
illustrated in Fig [6]), however the cross term vanishes. Nonetheless in the present case there 
is still a considerable synergy between the two broken symmetries. This shows up in the 
DOS as a pronounced dip at the chemical potential as illustrated in FigPD]. However, this 
effect is not confined to case of APdSC but is also present in the case of a uniform phase 
dSC coexisting with SDW as illustrated in Fig[9ti]. We have not found a simple explanation 
for this synergy. The quasiparticle bands in the vicinity of the Fermi energy have a complex 
dispersion for which we do not have a simple interpretation. Remarkably the form of the 
DOS near the Fermi energy is very similar for coexisting SDW and dSC with and without 
the array of vrDW the dSC This subtle difference in the DOS shows up as only a small 
difference in the ground state energy so that the energy cost of introducing vrDW is very 
small. 

VII. CONCLUSIONS 

The small energy difference that we find agrees with the earlier calculations reported by 
Himeda et air for coexisting SDW and APdSC^/^. These authors used a VMC method in 
which the strong coupling onsite constraint is exactly treated whereas here it is only ap- 
proximated through the Gutzwiller factors. This suggests that our failure to find a clear 
explanation for the stabilization of APdSC*/'' does not result from the Gutzwiller approxima- 
tion but may be because the t-J model omits some relevant physical effect. Alternatively the 
special cross term between SDW and APdSC order found by Agterberg and TsunetsugiJ^ 
which favors oblique wavevectors for the two periodicities may mean that our simple pattern 
with parallel arrays of AFDW and ttDW is not optimal, although on the surface it looks very 
plausible to simple stagger the two domain walls. After completing this paper, we learned 
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that a related work was posted by Chou et al.—. 
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FIG. 2: (Color online) The spatial (I) and wavevector (ky) dependence of the quasiparticle spectral 
weight Aj^ii.y{E) for the simplified model {Hs in Eq[TO]) with an isolated site-centered anti-phase 
domain wall in a d-wave SC. The parameters are the same as that used in Fig[l]. Panels (a) and 
(b) are for Ap = 0.02 (0.08), respectively. The energies E corresponds to the Andreev bound states 
(ABS) in the r.h.s. panels and the bulk SC gap in the l.h.s. panels as shown in Fig[T^,b]. In panel 
(a2) ABS extends away from the domain wall at site 28 into the bulk of the superconducting state 
due to |A| /Ep << 1, while in panel(b2) where |A| is much larger the ABS is much more confined 
in a small region around the domain wall. For the states close to the SC gap, small A leads to 
a more homogeneous state, while moderate A results in a great suppression of the state close to 
domain wall. 
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(a) A modulation with AP 
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FIG. 3: (Color online) (a) Schematic illustration of the modulation of the pairing parameter A for 
the simplified model (Hamiltonian Hg defined in Eq[lO]) for dSC state with periodic site-centered 
anti-phase domain walls with the shortest periodicity of Lx = 4. The anti-phase pattern of A 
is illustrated by the color scheme, (b) Local density of states (LDOS) with parameter values as 
doping concentration 5 = and Xp = ^p = 0.2. The domain walls are close to each other, they 
form bands with weak dispersion along kx but strong dispersion along ky parallel to the domain 
wall. At half filling, these bands display an anisotropic nodal structure as demonstated in Eq[T3] 
and by the low energy LDOS behavior. 
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(a.l) A modulation without AP 
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FIG. 4: (Color online) Schematic illustration of the modulations of the parameters pairing am- 
plitude A, hopping amplitude X) hole concentration S and antiferromagnetic moment m, for two 
states: SDW+CDW+dSC" [panels (a)] and SDW+CDW+APdSC* [panels (b)] (without and with 
site-centered anti-phase domain walls). The average doping is 1/8 and the periodicity L^ = 8. 
In panels (a/b 1-2) The amplitudes A and x ^I's denoted by the width of the bond, the spatial 
modulation of the staggered antiferromagnetic moment rrii is denoted by the arrows, the hole con- 
centration modulation is represented by the size of the dots. The anti-phasing of A in panel(b.l) 
is illustrated by the different color pattern at either side of the domain wall with cyan (magenta) 
for positive (negative) value. D-wave pairing symmetry is still preserved between two neighbor- 
ing domain walls. The anti-phase domain walls coincidence with the sites which have the largest 
sublattice antiferromagentic moment and smallest hole concentration. However, for the case with- 
out SDW, the domain walls locate at the sites with the largest hole concentratiouj^ Panels (a/b 
3) show the spatial hole density (red solid) and the AF moment (green dash) modulations. The 
site-centered anti-phase domain walls lead to an anisotropy of X) and an enhancement of the hole 
density and antiferromagnetic moment modulation. 
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FIG. 5: (Color online) (a) The energy (shown in Eq[5]) dependence of the two states 
SDW+CDW+dSC" and SDW+CDW+APdSC" (without and with site-centered anti-phase do- 
main walls) on the nnn hopping integral t' for isotropic and anisotropic nn hopping ratio tx/ty. 
The energy unit is to = 300mey. The nnn hopping integral t' does not but anisotropic t^(y) and 
Jx(y) do help to push the energy of SDW+CDW+APdSC'^ state (the solid and red symbol) closer 
to SDW+CDW+dSC* state (the open and blue symbol). Square (circle) symbols are for the values 
tx/ty = 1.00(0.85). Panels (b, c, d): the energy, charge and magnetization moment modulations 
of these two states with additional external potentials which are imposed to enhance the charge 
and magnetic modulations by shifting the local potential by +V (V > 0) up for the sites with 
the largest antiferromagnetic moment and —V down for the sites with zero antiferromagnetic mo- 
ment. A substantial anisotropy t^/ty = 0.85 is used. The diamond (triangle) symbols are for 
t' = — 0.3(— 0.1). The larger antiferromagnetic moment and hole concentration modulations drive 
the energy difference smaller between the two states. 
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(a. 1) A modulation without AP 



(b.l) A modulation with AP 
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FIG. 6: (Color online) Schematic illustration of the modulations of the parameters pairing am- 
plitude A, hopping amplitude X) hole concentration b and antiferromagnetic moment m, for two 
states: SDW+CDW+dSC'' [panels (a)] and SDW+CDW+APdSC'' [panels (b)] (without or with 
bond-centered anti-phase domain walls). The average doping is 1/8 and the periodicity Lx = 10. 
As shown in panel (bl) the anti-phase modulation of the pairing A is bond-centered with the 
domain wall located at the bonds connecting two nn sites with maximum stagger antiferromag- 
netic moment ib|r?T,| along x direction. The energy difference between the two states with and 
without bond-centered domain wall is even smaller than the case for site-centered domain wall 
with anisotropic nn hopping t^jty = 0.85. The modulation magnitude of the hole density and 
antiferromagnetic moment in these two states are close to each other. 
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FIG. 7: (Color online) dSC state with site-centered anti-phase domain wall but without antiferro- 
magnetism (doping concentration 6 = 1/8, t' = 0.0, V = 0, and a supercell Lx = 8, the energy unit 
is to = SOOmeV). The pattern for the pairing amplitude A is similar to the case shown in Fig [4]. 
(a) Fermi surface in the reduced Brillouin zone, (b) Quasiparticle (QP) DOS (blue) and average 
DOS (red). The two peaks A and B at negative energies are a consequence of the flat dispersion 
along ky direction [shown in panel (c)] formed by the propagating of Andreev bound state along 
y direction, (d) Local DOS (LDOS), near to the Fermi level the largest portion of the density of 
states is at the center of the domain wall. 
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FIG. 8: (Color online) SDW state (without dSC A = 0) with a periodicity of L^ = 8 and an 
average doping concentration 6 = 1/8, t' = —0.10 (the energy is in unit of to = SOOmeV). The 
antiferromagnetic sublattice moment pattern is the same as that shown in Fig[4]. (a) Fermi surface 
in reduced Brillouin Zone, (b) Local density of states (LDOS). Note that this SDW state is a 
metallic state. 
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FIG. 9: (Color online) SDW+CDW+dSC" and SDW+CDW+APdSC" states (without or with 
site-centered anti-phase domain wall). The upper figures (al, bl) show the local density of states 
(LDOS) at the three inequivalent sites with max \m\, zero \m\, and the middle site, (doping 
6 = 1/8, t' = —0.1, V = 0, and isotropic tx = ty). The energy is in unit of to = SOOmeF. The 
lower figures (a2, b2) show the average DOS and quasiparticle (QP) DOS. In order to facilitate the 
comparison between the states with and without the domain wall in panel (b2) the cyan curve is the 
QP DOS for the SDW+CDW+dSC* state, replotted from panel (a2). A small gap opens at zero 
energy. However, a substantial part of the DOS located at lower energy is pushed to closer to the 
Fermi level. This may be the reason that the opening of a gap in the SCW+CDW+APdSC* state 
does not lead to a lower energy relative to the state without the domain wall. Note that a small 
broadening 5 = 0.004ty is used to smooth the curve. The nodal behavior in SDW+CDW+dSC* 
state is not a general phenomenon. For larger t', anisotropic t^iy), or external additional potential 
this nodal structure may disappear. Also for other cases, (e.g. i' = 0, t^; = Ij ^ = 0) no gap opens 
in SDW+CDW+APdSC" state. 
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